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Abstract 

The Hu-Paz-Zhang equation is a master equation for an oscillator coupled 
to a linear passive bath. It is exact within the assumption that the oscillator 
and bath are initially uncoupled Here an exact general solution is obtained 
in the form of an expression for the Wigner function at time t in terms of the 
initial Wigner function. The result is applied to the motion of a Gaussian 
wave packet and to that of a pair of such wave packets. A serious divergence 
arising from the assumption of an initially uncoupled state is found to be due 
to the zero-point oscillations of the bath and not removed in a cutoff model. 
As a consequence, worthwhile results for the equation can only be obtained 
in the high temperature limit, where zero-point oscillations are neglected. In 
that limit closed form expressions for wave packet spreading and attenuation 
of coherence are obtained. These results agree within a numerical factor with 
those appearing in the literature, which apply for the case of a particle at 
zero temperature that is suddenly coupled to a bath at high temperature. On 
the other hand very different results are obtained for the physically consistent 
case in which the initial particle temperature is arranged to coincide with that 
of the bath. 
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I. INTRODUCTION 



The Hu-Paz-Zhang equation is a master equation with time-dependent coefficients for a 
harmonic oscillator interacting with a linear passive heat bath of oscillators. The equation is 
exact and general within the assumption that in the initial state the bath is in equilibrium 
and not coupled to the oscillator. It was ffist derived in generality using path integral 
methods by Hu, Paz and Zhang [1], although equivalent equations had been obtained earlier 
for the case of an Ohmic bath [2,3]. See also [4]. Later a derivation in the form of an equation 
for the Wigner function was given by Halliwell and Yu [5] , who corrected a misprint in the 
earlier publication. Using the notation of these last authors, the equation has the form 

dW 1 dW X dW 

at m oq op 

dvW d'^W d'^W 

+2r(t)^ + hmmhit)^ + ^r(t)/w^, (1.1) 

where ^'^{t), 2r(t), h{t), and f(t) arc time-dependent parameters for which one has explicit 
expressions (see equations (2.19) and (3.8) below). 

The integration of this equation, with its time-dependent coefficients, appears to be a 
formidable problem. Indeed, earlier discussions have generally used numerical methods. Our 
purpose here is, ffist of all, to present an exact general solution of this equation. This solution 
can be evaluated in explicit, closed form for many problems of interest. In particular, we 
exhibit the solution for two such problems: an initial state corresponding to a Gaussian 
minimum uncertainty wave packet and an initial state corresponding to a widely separated 
pair of such wave packets. We use these results to accomplish our second purpose, which 
is to critically examine the assumption of an uncoupled initial state. We find that there 
is a serious difficulty arising from this assumption: the zero-point oscillations of the bath 
give rise to a divergence that leads to an instantaneous spread of a wave packet to infinite 
width. In effect, the state instantaneously disappears! The result is that meaningful results 
can be obtained only in the high temperature limit, where one conventionally neglects the 
zero-point oscillations. Even in this limit, we find significant difficulties arising from the fact 
that translational invariance is broken. Nevertheless, we find for short times expressions for 
wave packet spreading and attenuation of coherence that are consistent with those found 
by earlier authors. On the other hand, by adjusting the initial temperature of the particle 
to be the same as that of the bath, we find in the high temperature Ohmic limit results 
consistent with exact calculations which take into account entanglement at all times [6,7]. 

The plan of this paper is as follows. The basis for our discussion is the quantum Langevin 
equation, so we begin in Section II with a description of that equation and its solution, first 
for the stationary case, then for the initial value case, and finally for the form local in time. 
We next in a short Section III we give a simple derivation of the exact master equation 
(1.1). Then in Section IV we derive our general solution. The key result, given in (4.15), 
is an explicit expression for the Wigner function at time t in terms of the initial Wigner 
function. A particularly useful result, given in (4.24), is an expression for the probability 
distribution at time t. In Section V we first evaluate this expression to find the probability 
distribution corresponding to an initial Gaussian wave packet. There we find the divergence 
mentioned above, which in Appendix A is shown to be present even in a model with a 
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high frequency cutoff. Also in section V we consider the motion of a pair of Gaussian wave 
packets (Schrodinger "cat" state) and obtain an exphcit expression for the attenuation of 
coherence. For the case of a particle at temperature zero suddenly coupled to a bath at high 
temperature T, this leads to an expression for the decoherence time equivalent with that 
appearing in many places in the literature. But for a consistent initial state in which the 
temperature of the particle is adjusted to coincide with that of the bath, a very different 
expression for the decoherence time is obtained, an expression corresponding to decoherence 
without dissipation [7]. Finally, in Section VI we summarize our results and make some 
concluding remarks. 



II. THE LANGEVIN EQUATION 

A. Stationary process 

The Langevin equation is a Heisenberg equation of motion for x{t), the dynamical vari- 
able corresponding to the coordinate of a Brownian particle in equilibrium with a linear 
passive heat bath. For the case of a particle in an external oscillator potential, this equation 
for the stationary process has the well known form [8], 

mx+ f dt'ii{t-t')x{t') + Kx^F{t), (2.1) 

where /x(t) is the memory function, K is the oscillator force constant and F{t) is a fluctuating 
operator force with mean {F(t)) = 0, and whose correlation and commutator are given by 



1 1 f°° huj 

- (F(t')F(t) + F(t)F(t')) =- / dujRe{jl(u + iO'^)}nw coth-— cos u;(t-t'), 

2 TT JO 2kT 

2?7, roo 

[F(t),F(t')]^ — dujRe{jl(LJ + iO+)}uj sin uj(t-t'). (2.2) 

ITT Jo 

Here il{z) is the Fourier transform of the memory function, 

roo 

fl(z) = / dtn{t)e'^\ (2.3) 
Jo 

As a consequence of the second law of thermodynamics, ji{z) must be what is called a positive 
real function: analytic and with real part positive in the upper half plane. In particular, 
jl{z) can be represented in terms of the real part of its boundary value on the real axis 
through the Stieltjes inversion theorem 

fliz) = -icz + — / duj^^^ 2.4 

TT Jo z"^ — o;^ 

where c is a positive constant. 

The solution of the Langevin equation (2.1) can be written 

x^{t)= f dt'G{t-t')F{t'), (2.5) 

J— CO 
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where G{t), the Green function, is given by 

1 

G{t) = — duja{uj + iQ+)e-'''\ (2.6) 
with a{z) the famihar response function 

«(^) = 2 ^-(^^T^ - (2-7) 

— mz^ — zzii[z) + K 

Here we have introduced a subscript s to emphasize that Xs{t) is a stationary operator- 
process, in the sense that correlations, probabihty distributions, etc. for this dynamical 
variable are invariant under time-translation {t ^ t + to). In particular, the correlation. 



{xs{t)xs{t') + Xs{t')xs{t)) = — d(jjlm{a{u! + iO"*")} coth — — cosa;(t — t'), 

TT Jo 2kT ^2.8) 



is a function only of the time-difference t — t'. In addition, for the free particle, where K = 0, 
the process is invariant under space translation {x —>■ x + a). 



B. Langevin equation for the initial value problem 

The description of the system given by the Langevin equation can be realized by a bath 
of harmonic oscillators. Perhaps the simplest such system, and the one we use as the basis of 
our discussion of the Hu-Paz-Zhang equation and its solution, is the independent oscillator 
model, for which the Hamiltonian is [9,8] 

^-|^+i^-^+E{|-+i™.-ite-n. (2.9) 

Writing the equations of motion and then eliminating the bath variables in terms of their 
initial values, one obtains the Langevin equation for the oscillator with given initial values 
[10], 

mx + /* di!ii{t - t')x(t') + Kx^ -fi(t)x(0) + F(t), (2.10) 
Jo 

where the memory function is given by 

H{t) = ^rrijUj'^j cos ujjte{t), (2.11) 

3 

while the random force is given in terms of the initial bath variables by 

F{t) = ^{gj(0)mja;| cosa;^^ + Pj{^)<^j sin oujt}. (2-12) 
j 

To express the solution of this equation, we first note that the Green function (2.6) van- 
ishes for negative times and for positive times is the solution of the homogeneous equation, 

mG + f dt'ii(t - t')G(t') +KG^O, (2. 13) 

^0 
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with the initial conditions 



G(0) = 0, G(0) = — . (2.14) 

m 

With this, we can show that the general solution of the initial value Langevin equation (2.10) 
is given by 

x{t) = mG{t)x{0) + mG{t)x{0) + X{t) 

x{t) = mG{t)x{0) + mG{t)x{0) + X{t), (2.15) 
where we have introduced the fluctuating position operator, 

Xit)= f dt'Git-t')Fit'). (2.16) 

^0 

In our subsequent discussion we assume that at i = the system is in a state in which 
the oscillator is not coupled to the bath and that the bath is in equilibrium at temperature 
T. In particular this means that the initial coordinates of the oscillator are not correlated 
with those of the bath, i.e., (x(O)-F(t)) = {x{{))F{t)) = 0. On the other hand, with regard 
to the bath, the equilibrium is with respect to the bath Hamiltonian, i^bath, obtained by 
setting the oscillator variables x and p equal to zero in (2.9), 



^bath = E(^ + ^^.^M)- (2-17) 



"J 



With this we find {F{t)) = 0, and the correlation and commutator are the same as those for 
the stationary equation, given in (2.2). 

Typically, the memory function fi{t) falls to zero in a very short time r, called the 
relaxation time of the bath. For times long compared with this bath relaxation time, the 
extra term on the right hand side of (2.10) vanishes, but only for much longer times, times 
long compared with the oscillator decay time, will this equation become the stationary 
equation, with the lower limit on the integration taken to be — oo. To be more specific, we 
note from the general expression (2.10) for the Green function that, so long as the oscillator 
force constant K is not zero, the Green function will vanish exponentially for long times. 
This follows from the Tauberian theorem: the asymptotic behavior of a function depends 
upon the low frequency behavior of its Fourier transform. It follows that, for long times, 
the dependence upon the initial coordinates in (2.15) disappears and, from a comparison of 
(2.16) with the expression (2.5) for Xs{t), that X{t) becomes the solution of the stationary 
Langevin equation (2.1). 



C. Form local in time 

We want now to write the Langevin equation (2.10) in the form of an equation that is 
local in time with time-dependent coefficients. To get this form, we first invert the equations 
(2.15) to express the initial variables in terms of those at time t. We next form the time 
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derivative of the second of the equations (2.15) and then insert these expressions for the 
initial variables in the right hand side. We can write the result in the form 

x + 2T{t)x + n^{t)x = —F{t), (2.18) 

where we have introduced the quantities 

G{t)G{t) - G{t)G{t) dlog{G^ - GG) 



2r{t) = 



G^{t) - G{t)G{t) dt 



^ G\t)-G{t)G{t) .2 19) 

^ ' G^{t) - G{t)Git) ■ ^ ■ ^ 

This equation is the local form we seek. In obtaining this form we have used the fact 
that X{t) is the solution of the inhomogeneous equation (2.10) with the initial conditions 
X{0) — and ^(0) = 0. It follows that X{t) is also the solution of the local equation (2.18) 
with the same initial conditions. The fact that the symbols used for the quantities (2.19) 
also appear in the exact master equation (1.1) is not accidental, as we shall see in the next 
section, where we give a derivation of that equation. 



III. THE TIME-DEPENDENT MASTER EQUATION 

The strategy for deriving the exact master equation is to obtain expressions for the 
first and second moments, first from the equation (1.1) and then from the local Langevin 
equation (2.18). From a comparison, wc obtain explicit expressions for the time-dependent 
parameters in the Hu-Paz-Zhang equation (1.1). 

In forming the moments of the equation (1.1), we take note of the relations [11] 

, ih d , , ih d , 

xp^iQ + ^Q^W, px ^ {q- -—)W, 

Here on the left x and p are the position and momentum operators for the oscillator, while 
p is the density matrix operator. On the right, q and p are the c-number variables of the 
Wigner function W{q,p;t). Thus, for example, 

fOO POO jp. f) 

(x) = Tr{xp} = dq dp{q + -t-)W. (3.2) 

J-oo J-oo 2 op 

In this way, forming the first moments of equation (1.1), we find 

(x) = ^ (P) , (P) = -2r(t) (p) - mQ\t) (x) . (3.3) 
Eliminating (p), we find 

{x) + 2T{t) {x) + n\t) {x) = 0. (3.4) 
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If now, we form the mean of the local Langevin equation (2.18). using the fact that {X{t)) = 
0, wc get exactly the same equation, but with now the quantities 2r(t) and fl'^{t) given by 
the expressions (2.19). 

Next, forming the second moments of the equation (1.1), we find 

d{x^) 1 , , 
dt m 

'^^''P+P'') ^-/p')- 2mn' (x') - 2r {xp + px) + 2hT /, 
dt m \ / \ ' 

^ ^ ^ - -mfi^ {xp + px) - 4r + 2hmTh. (3.5) 



dt 

On the other hand, using the local Langevin equation (2.18), we find 

d{x'') 



dt 

d {xx + xx) 
di 



~ {xx + xx) , 

= 2 ^i^^ + {xx + xx) 

= 2 - 2Q^(i) (x^) - 2r{t) {xx + xx 

+-{x{t)F{t)+F{t)xit)), 
m 



d{x^) 
dt 



— {xx + xx) 

= -Q'^{t) {xx + xx) - 4r(t) 

+- {x{t)F{t) + F{t)x{t)) . (3.6) 
m 

Now, in the right hand side of the second of these equations we use the fact x{t) — X{t) is 
not correlated with F{t) This should be clear since, as we see from (2.15), this combination 
depends only upon the initial coordinates of the oscillator. Therefore we can replace x{t) 
with X{t) in the correlations of x with F. Using the same argument in the last of these 
equations, we see that after a little rearrangement we can write 



d{x^) 
dt 

d {xx + xx) 
di 



= {xx + xx) , 

= 2 (^x^) - 2Q\t) (x^) - 2V{t) {xx + xx) 

+-{X{t)F{t) + F{t)X{t)), 
m 

= -VL^it) {xx + xx) - AV{t) (x 

+-{xit)F{t) + Fit)X{t)). (3.7) 

We now compare these equations with the equations (3.5) obtained from the equation 
J^l.l). In doing so we must interpret p = mx. We see then that we can identify 



d{x'') 
dt 
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2hr{t)f{t) = {X{t)F{t) + F{t)X{t)), 
2hr(t)h(t) = (^X{t)F(t) + F{t)X{t)) . 



(3.8) 



This completes the derivation of the exact master equation (1.1), with exphcit expressions 
for the time-dependent coefficients. 



IV. GENERAL SOLUTION OF THE TIME-DEPENDENT MASTER EQUATION 

The task of solving the equation (1.1), with its time-dependent coefficients given by the 
complicated expressions (2.19) and (3.8), appears formidable. Indeed, if one were presented 
with the equation with no idea of the origin of the coefficients its solution would be very 
difficult. But we have in (2.15) an expUcit solution of Langevin equation describing the 
underlying motion. This will allow us to construct the general solution of the equation. 

To begin, we remind ourselves (1.1) is an equation for the reduced density matrix, given 
by the partial trace over the bath coordinates. That is, 

W{q,p;t) = j dqj rfpWsystem(g,p; q, p; t)- (4.1) 

Here VFsystem is the Wigner function for the system of oscillator and bath, with q = 
{(liyQ2,-'') ^-iid p = {pi,P2,-'') the bath coordinates and momenta. Now, the system 
is one of coupled oscillators and for such a system the solution of the equation of motion 
is formally identical to that for the corresponding classical system. That is, the Wigner 
function for the system at time t is related to that at time t = through the relation 

VVsystem(g,p;q,p;i) = W^system(g(0),p(0);q(0),p(0);0), (4.2) 

where g'(0),p(0); q(0), p(0) are the initial values for which the solution of the equations of 
motion is such that q{t) — q, p{t) = p, q(t) = q, p(t) = p. Finally, we remind ourselves 
that for the Hu-Paz- Zhang equation the initial state is a product state corresponding to a 
Wigner function of the form [1 1] 

Wsystem{q,p; q, p; 0) = W{q,p; 0)l[wj{qj,pj). (4.3) 

j 

Here, on the right W{q,p; 0) is the initial Wigner function for the oscillator and the product 
is the Wigner function for the bath, in which Wj{qj,pj) is the Wigner function for a single 
oscillator of mass rrij and frequency cuj, 

1 r p"^- + Tfi^-uj^q^ \ 

Combining these results, we see that the reduced density matrix at time t is given by 

W{q,p;t) ^JdciJ dpW{q{0),p{Oy,0)l[wj{qj{0),pj{0)). (4.5) 
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This reduced density matrix is the solution that wc seek. So far, however, aU we have done 
is to carefully indicate the definition of this quantity, we must now carry out the indicated 
operations to obtain an explicit expression. 

As a first step we transform the integration to the initial bath coordinates, holding q and 
p fixed. Under this transformation, 

^q^P = /^(^'^'^'P) dq(0)dp(0), (4.6) 

d{q,p; q(0),p(0)) 

where the factor is the Jacobean of the transformation, for which we have used the notation 
of Landau and Lifshitz [12]. But, 

%,P;q,p) _ d{q,p; q,p) (9(g(0), p(0); q(0), p(0)) 



d{q,p; q(0), p(0)) 9(g(0), p(0); q(0), p(0)) d{q,p; q(0), p(0)) 
_ %(0),p(0);q(0),p(0)) 
d{q,p; q(0),p(0)) 

^ / %,p;q(0),p(0)) y 
V%(0),p(0);q(0),p(0))y' 

_ ( dq dp dq dp \ 

\dq{0)dp{0) dp{0)dq{0) J ' ^ ' ^ 

Here, we use in the first line the fact that the Jacobean of two successive transformation 
is the product of the Jacobians, in the second line the fact that the motion of the system 
corresponds to a canonical transformation for which the Jacobean is unity, in the third line 
the fact that the Jacobean of the inverse transformation is the reciprocal of that of the direct 
transformation, and finally in the last line the definition of the Jacobean as the determinant 
of the matrix of partial derivatives. Now, to evaluate this Jacobean we use the solution 
(2.15), which we write in the form 

q = q{t) = mG{t)q{0) + G{t)p{0) + X{t), 

p = p{t) = m^G{t)q{Q) + mG(t)p(0) + mX(t). (4.8) 

Here we recall, from the definitions (2.16) of X{t) and the expression (2.12) for F{t) that 
X{t) depends only on the initial coordinates of the bath, which are held fixed in forming 
the partial derivatives in the last line of (4.7). Therefore, we see that 

%,p;q,p) 1 ^4_g^ 



%,p;q(0),p(0)) m'^iG^-GG) 
and, using (4.6) we can write (4.5) in the form 

>^fa.P;^)^ ™-'"^'> (4.10, 

where the brackets represent the average over the initial equilibrium distribution of the bath. 
Again, we remind ourselves that in the integrand g(0) and p(0) are obtained by inverting 
the equations (4.8). That is. 
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_ mG{q-X)-G{p-mX) 
' m?{G''-GG) 
^ -m'Gjq -X)+ mGjp - mX) 
^ ' m?{G-' -GG) ' ^ ' ' 

Since X is linear in the initial bath variables, its average has the Gaussian property: averages 
of all moments can be expressed in terms of those of the second moment. 

We can carry out this average if we introduce the Fourier transform of the initial Wigner 
function, writing 

= (2^ loo'^^ y_^rfPl^(g,P;0)e'(^«+«^)/'^ (4.12) 
Inserting this in (4.10), we can write 

1 pCO poo 

WiQ,P;t) = ,o . ^2fn2 nn^ / / dPW{Q,P-Q) (exp{z[Pg(0) + gp(0)]/n}) . 

[27rnm)^{G^ — GG) J-oo J-oo ^^-^ 

The form of this result can be made a bit simpler if we introduce a transformation to 
variables r and s, defined by 

Q ^ mGr + Gs, P = rn^Gr + mGs. (4. 14) 

We then find that Pg(0) + Qp{Q) = r{p - mX) + s{q - X) and dQdP - m'^{G'^ - GG)drds, 
so that (4.13) becomes 

2 poo poo ^ . ... 

W{q, p; t) = - — — / dr dsW{mGr + Gs, m^Gr + mGs; 0) 

(^TT/lj J —oo J ~oo 

where we have used the Gaussian property to write 

/^-i{mXr+Xs)/!i\ ^ ^-^{m^X^y+m{xX+XX)rs+{x^)s^)_ ^^_-^g^ 



This is the form of the solution that is perhaps most useful. In its evaluation, the Green 
function G{t) is given by (2.6), while X{t) is given by (2.16) and its correlations are evaluated 
using (2.2). 

While we find (4.15) to be the most useful form of the solution, one can insert the inverse 
of the Fourier transform (4.12) to express the solution in the form of a transition operator 
acting on the initial Wigner function, 

/oo poo 
dq' / dp'P(?,p;g',p';i)iy(g',p';0). (4.17) 
-co J—oo 

Here P{q,p;q',p';t), called the transition probability, can be written 

P{q,p;q',p';t) = — -L==exp{-lR- A-i.R}, (4.18) 
27rvdetA 2 
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where we have used a dyadic notation with 



' f (jf A' + XX) (A'2) 

H«)^(::l^i:j>). (4.19) 

Here, in R, the quantities {q(t)) and {p{t)) correspond to the mean of the initial value 
solution (2.15) with initial values q' and p'. That is, 

{q{t)) = mG{t)q' + G{t)p' , 

{p{t)) ^m^G{t)q' + mG{t)p'. (4.20) 

It is of interest that this expression (4.18) for the transition probabihty is formally the 
same as that for the classical Kramers equation [13]. The difference is that the Green 
function and the mean square of the fluctuating position and velocity operators here are for 
a quantum oscillator interacting with an arbitrary heat bath, while in the classical solution 
of the Kramers equation they are for a classical oscillator interacting with an Ohmic bath. 
Another significant difference is that the solution (4.20) is that of the mean of the initial value 
Langevin equation (2.10), with the term — //(i)x(O) on the right hand side. The classical 
solution of the Kramers equation corresponds to dropping this term, since it is describing 
the classical stationary process. 

As a first illustration of the utility of the form (4.15), we show how the equilibrium 
solution arises for long times. First we recall that, so long as the oscillator force constant K 
is not zero, the Green function will vanish as i — > oo. Next, we recall that, again as i — > oo, 
X{t) Xs{t), the solution (2.5) of the stationary Langevin equation (2.1). Thus, we see 
almost by inspection that 



where (Xg) and (ig) are the equilibrium values of the mean square position and velocity, 
given by the fluctuation-dissipation theorem. 



xz) —— (i6(;Im{Q;(a; + iO"*")} coth 
' 71 Jo 



71 Jo ^ ^ 2kT' 

- / o?u;u;2lm{a(a; + i0+)}coth— — . (4.22) 
TT Jo 2kT 



This result is perhaps more familiar in the weak coupling limit, where Im{Q;(a; + iO"^)} — > 
Tr5{uj-uJo)/2mwo with cuq = sjK/m. Then {xfj = ul {xl) = ^ coth ^ and (4.21) becomes 
the familiar equilibrium form of the Wigner function for the uncoupled oscillator [11]. 

Finally, we remark that the interest is most often in the probability density at time t, 
given by 

/oc 
dpW{x,p;t). (4.23) 
-oo 
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Using the solution (4.15), the integral over p gives a 5-function in r. With this we can 
perform the r integration to obtain the result 

P{x; t) = ^ dsW{Gs, mGs- 0) exp{-^ (x^) + i-s}. (4.24) 

In the next Section, we evaluate this probability density for some problems of interest. 



V. EXAMPLES 

In this section we consider the evaluation of the general result (4.24) for the case of a free 
particle {K — 0) interacting with an Ohmic bath. We have chosen the Ohmic model since 
it is that used most extensively, almost universally, in discussions of dissipative systems. 
Newtonian drag is Ohmic, as is the Stokes force and, of course, classical Brownian motion. 
In quantum electrodynamics, the Weisskopf-Wigner approximation is an Ohmic model. In 
addition, our results take their simplest form for that model. In Appendix A we present 
selected results for the more general single relaxation time model. 

The examples are intended, first of all, to illustrate the power and utility of our exact 
solution. They arc chosen since they appear in a truly large number of recent papers, 
where approximate methods were used. They are also the examples discussed, by numerical 
methods, in the original Hu-Paz-Zhang paper as well as in the earlier papers we have cited. 



A. Preliminary formulas 

For the Ohmic bath the memory function has the form 

li{t) = 2C5(t), (5.1) 

where ( is the Newtonian friction constant (the factor 2 is because the integral in the 
Langevin equation is over only half the delta-function). In this case the equation (2.10) is 
already in local form. With the form (5.1) for the memory function and with K = 0, the 
response function (2.7) takes the simple form 

a{z) = (5.2) 



The Green function (2.6) is then 



1 _ p-Ct/m 

G{t) = . (5.3) 



If we form the quantities (2.19) with this Green function, we find 

2r{t) = ^, n\t) = —s{t). (5.4) 

m m 

For any finite time these expressions follow trivially. The delta function is not so easy to see, 
although it should be obvious from the form (2.10) of the Langevin equation. In Appendix 
A, this result is derived explicitly in the Ohmic limit of the single relaxation time model. 
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The only other quantity that we need to evaluate the general result is {X'^{t)). the 
mean square of the fluctuating position operator. In this connection, it is a simple matter, 
comparing the stationary solution (2.5) with the definition (2.16) of X{t), to obtain the 
general relation 

X{t) = Xs{t) - x,{0) + f dt'{G{-t') - G{t - t')}F{t'). (5.5) 

— oo 

In discussing this operator, in particular the mean of its square, we shall make use of the 
mean square displacement for the stationary process, 

s{t) = {[x,{t) - xM?) 

'2fl fOO fyjj 

= — / du!lm{a(u! + iO'^)} coth — — (l — coscut). (5.6) 
TT JO 2kT 

The discussion in the previous paragraph has been general, applying to an oscillator 
interacting with an arbitrary heat bath. We now specialize to the current case of a free 
particle interacting with an Ohmic bath. The Green function is then given by (5.3), from 
which we see that G{—t') — G{t — t') — —mG{t)G{—t') and (5.5) becomes 

X{t) = Xs{t) - Xs(0) - mG{t)xs{0). (5.7) 

Forming the mean square, we can write 

[X\t)) = s{t) - mG{t)s{t) + ^m^G\t)s{0). (5.8) 
In the same way we find 

X^t)) = ^[1 + m^G^]s{0) - mG{t)s{t). (5.9) 

For this Ohmic case, using the form (5.2) for the response function in the expression (5.6) 
for the mean square displacement, we find 

s(t) = ^ /"^^,^,,, (l-cosa;t). (5.10) 

1. High temperature limit 

In the high temperature limit {kT ^ h(/m) we replace the hyperbolic cotangent in 
(5.10) by the reciprocal of its argument (thus, neglecting the zero-point oscillations) The 
result takes the form [14] 

2kT , l-e-^*/"*, 
s{t)^^{t-m ). (5.11) 

With this in the expression (5.8) for the mean square of the fiuctuating displacement, we 
see that, still in the high temperature limit. 



X\t)) ^ —{2t - 2m '-). (5.12) 
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2. Zero temperature 



At zero temperature, we replace the hyperbolic cotangent in (5.10) by unity. The result 
can be written 



s{t) = ^/(^), (5.13) 



where [14] 



I{x) = / dy (1 - cosy) 

Jo y{y^ + x^) 

= logx + 7 - ^[e-^Ei(x) + e^Ei(-x)]. (5.14) 
Here 7 = 0.577215665 is Euler's constant. Note the expansions [15], for small x, 

1 00 n 00 / \n 

I{x) = -{logx + 7)(coshx - 1) - -[e-- ^ ^ + ^ L_^], (5.15) 

2 ;S ri\n 

and, asymptotically, for large x, 

I{x) ~ log X + 7 - ^ - ^ - ^ - • • • . (5.16) 

Here we see that there is a serious concern: for small x the second derivative I"{x) = — log a; 
and therefore the term involving s(0) in the expression (5.8) for {X'^{t)) is logarithmically 
divergent. This divergence persists for long times, where 

X'{t))^^\ogCt-^\ogO\ (5.17) 

in which the neglected quantity is of the order of a finite constant. This divergence for the 
Ohmic case has, of course, been noted by earlier authors [3,1], but it docs not seem to be 
known that this divergence persists in a model with a high frequency cutoff. In Appendix 
A we show this explicitly for the single relaxation time model. 



B. Gaussian wave packet 

To begin, we note that in evaluating the solution we make repeated use of the standard 
Gaussian integral: 

' 2 



roo i /27r 6 

/ dxexp{ — ax'^ + ibx} = \ — exp{ }. (5.18) 

J-00 2 V a 2a 



We consider an initial state corresponding to a Gaussian wave packet of the form. 
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This is a so-called minimum uncertainty wave packet, centered at {x{0)) — xo and 
with width (Aa:^(0)) = ([a;(0) — {x{0))]'^) = cr^. The corresponding mean momentum is 
(p(0)) = m {x{0)) = and the momentum width (corresponding to minimum uncertainty) 
is (A]9^(0)) = — (2^(0))]^) = The Wigner function corresponding to this state 

is 

1 /"^ ?/ 11 

W{q,p; 0) = ^ due^-^/'^ijiq - -, 0)^*(g + -, 0) 

1 7 (g-^o)' 2aV, 

Its Fourier transform is 

W{Q,P;0)^ dq dpe-'^^''+^P^/''W{q,p;0) 
J— 00 J —00 



L 



—00 
00 



(i,e-"''/Vta-|,Om<;+|,0) 



= exp{-«;-:^-.^|^}. (5.21) 
Putting this in the expression (4.24) for the probability density at time t, we find 

P(x; t) = , exp{- ^ ^.^ }, 5.22 

^ ^ ^27r(Ax2(t)) ^ 2(Aa;2(i)) ^' ^ ^ 

where (Aa;^(t)) is the variance of the position, given by 

{Ax\t)) = m'G\t)a' + + (^^(i)) . (5.23) 

This is a general result, valid for any heat bath. For the Ohmic case the Green function 
is given in (5.3) and {X'^{t)) is given by (5.8). For the single relaxation time model, the 
corresponding results are given in (A4) and (A9) 

The first thing that should strike us in this result is that the variance is in fact infinite, 
due to the divergence found in the expression (5.8) for {X'^it)) at zero temperature. This 
divergence arises from the zero-point oscillations and is therefore always present, although 
conventionally one neglects the zero-point oscillations in the high temperature limit. This 
is a serious difficulty, since it tells us that an initial wave packet spreads instantly to infinite 
width. In effect, the state vanishes! We emphasize that the problem here is not with the 
Hu-Paz-Zhang equation or its solution, which are perfectly correct, but with the assump- 
tion of an uncorrelated initial state. Unease with this assumption has been expressed by 
many authors (we note in particular the remarks in the concluding Discussion section of the 
Hu-Paz-Zhang paper [1]) but it does not seem to have been realized how serious are its con- 
sequences. Indeed, the only meaningful results for the equation are for the high temperature 
limit, and we confine our discussion in the following to that limit. 

In the high temperature limit, using the expression (5.3) for the Green function in the 
Ohmic case, the probability density (5.22) becomes 
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where, evaluating the expression (5.8) for {X^{t)) with the high-temperature expression 
(5.11) for the mean square displacement, the variance is given by 



/a 2. n\ 2 2Ct/m fv'il - e'W^f 

+— (2t-2m '-). (5.25) 

A difficulty with this result is that the center of the wave packet drifts to the origin. Since 
for a free particle the origin cannot be a special point, we see from this that the translational 
invariance of the problem is broken by the assumption that the initial state corresponds to 
an uncoupled system. Indeed, the system Hamiltonian (2.9) is invariant under simultaneous 
translation of the particle and bath coordinates (a; — > a; + qj — > qj + d) and the time- 
dependent master equation (1.1) correctly describes the dynamics of the system with regard 
to this Hamiltonian. But the bath Hamiltonian (2.17) does not possess this invariance and 
the initial state is one in which the bath is in equilibrium with respect to this Hamiltonian. 
Another way to see this effect is to note that for a free particle interacting with an Ohmic 
bath the mean motion described (1.1) satisfies the equation 

m(x) + C(i) = -2C5(t)(x(0)). (5.26) 

That is, the particle receives an initial impulse — C(a;(0)), directed toward the origin and 
with a magnitude such that in the course of time the particle arrives at the origin. Another 
difficulty, which in fact has the same origin as the first, is that in the expression (5.25) for 
the variance the first term, which corresponds to the initial variance, decays in time. But 
the initial variance should persist and not decay. 

A further difficulty is seen if we look at the variance (5.25) for times short compared 
with the Ohmic decay time , 

(Ax^(t))-a^ + ^, a/m«l. (5.27) 

But this is exactly the formula for the RMS width of a wave packet as obtained from 
elementary quantum mechanics [16]. It therefore corresponds to temperature zero and the 
thermal spreading one should expect at high temperature is absent [17,6,7]. Now this last is 
not a difficulty with the equation or our solution, but with the initial state we have chosen, 
which corresponds to a particle at temperature zero. Instead we should choose a state 
corresponding to a wave packet at temperature T, obtained by averaging the initial Wigner 
function (5.20) over a thermal distribution of initial velocities. 

To accomplish this, we note first that the initial state (5.19) corresponds to a particle 
at rest (i.e., at T = 0). To obtain the state corresponding to a particle with velocity v, we 
simply multiply with a factor e*™"^/'*. With this, we see that the Fourier transform of the 
corresponding Wigner function is obtained by multiplying (5.21) by e"*""'*^/'^. The thermal 
average of this factor is 

m f°° mv'^ mQ 

e^Pi-TTI^ - ^^^) = exp{-;rYi}, (5.28) 



2nkT J-oo 2kT h ' ^ ' 2A 
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where A is the thermal de Broghe wavelength, 

A=^^. (5.29) 

Therefore, we see that for a particle at temperature T the Fourier transform of the initial 
Wigner function is obtained by multiplying the corresponding function at T = by the 
factor (5.28). Multiplying (5.21) by this factor we obtain 

P; 0) = exp{-(^ + ^)Q^ " W " '^>- (^■^•" 

where we indicate that this corresponds to temperature T by the subscript. Using this in 
(4.24) we find for the probability distribution at temperature T the result 

where we have introduced 

A.%=(A.^) + ^(l-e-^^/-)^ 

= <T^e-^*/- + ^^^^'"^ + ^X\t))^ , (5.32) 

in which 

X'it))^ = ^—{2t - 2m^ ). (5.33) 

Note that now the variance for very short times includes the thermal spreading [6,7], 

Mt-'^' + Z^ + — C^/m«l. (5.34) 

However, the long time drift of the wave packet center and the shrinking of the initial 
variance remain. 



C. Pair of Gaussian wave packets 

Wc consider now an initial state corresponding to two separated Gaussian wave packets. 
The corresponding wave function has the form 

^(^, 0) = (exp{-^^^} + exp{-^^f|^}), 

where d is the separation and a is the width of each packet. Using (5.21) we see that the 
Fourier transform of the initial Wigner function is given by 

W(Q,P;0) = _^exp{-^ - _}(oos- + e'i;^ cosh 



e 8<T^ 



(5.36) 
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Putting this in the expression (4.24) for the probabihty density at time t, we find 

1 / {x - tC^'^''"')'^ ( + -C//™)2 



2 



^2 (1 _ g-Ct/m)2;^2^2 



Here the first two terms in the parentheses correspond to a pair of Gaussian wave packets 
of the form (5.22), initially centered at x = ±d/2 and drifting toward the origin. The third 
term, that with the cosine, is an interference term. The attenuation factor a{t) is the ratio of 
the coefficient of the cosine term divide by twice the geometric mean of the first two terms. 
We find 

"(t^^'M- Jf^^y^ -ei (6.38) 

where at high temperature {X^{t)) is given in (5.12) and {Ax'^{t)) in (5.25). 

For times long compared with the Ohmic decay time, t ^ rn/(, we see that (X^(t)) 
and (Ax^) become asymptotically equal, growing linearly with t. For such long times the 
attenuation coefficient (5.38) therefore approaches a very small constant, a(t) ~ exp{ — ^}. 
For times short compared with the Ohmic decay time, t -C Tn/(, we see that {X'^(t)) = ^^^^ 
while {Ax'^{t)) = a"^ + and therefore 

If we suppose that the slit width is negligibly small, we find a{t) = exp{— t/r^i} where 
Ta = This, except for a factor of 6 is exactly the decoherence time that often appears 

in the hterature [18]. But, as we have seen above, this result corresponds to a particle in 
an initial state that is effectively at temperature zero, which is suddenly coupled to a heat 
bath at high temperature. The result is therefore unphysical in the sense that the initial 
state does not correspond to that envisioned when we speak of a system at temperature T . 

Now, just as in our discussion of the spreading of a single Gaussian wave packet, this 
difficulty can be repaired by choosing the initial state of the particle to be the same as that of 
the heat bath. The prescription for accomplishing this is very simple: multiply W{Q, P; 0), 
given in (5.36), by the factor given in (5.28). The result is to replace (5.36) by 

_ + ^)n2 _ ^Lil 
8a^ 2~X^'^ 2^2 



W^iQ, P; 0) = ^ exp{-(— + ^)Q^ - ^} 



x(cos— + e 8^ cosh ^|-^). (5.40) 



1 + e Sa 

Pd _^ , Qd 
— - + 6 8^ cosh — 
2h Aa 

Putting this in the expression (4.24) for the probability density at time we find, in place 
of (5.37) 

P-rix; t) = , exp{--^ ^— ^ — -\ + exp{--^ ^— r — —\ 

^ ' 2(l + e-<^V2a^)y2^r(A^ V 2{Ax^)^ ^ 2{Ax^)^ ^ 

x^ {1 - e-W'^fh^d^ GMx \ , , 

^"P^- 2(A^ - 8^ + 32CVMAX2), ^ Aa- {Ax-), ) ' ^'-"'^ 
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where {Ax'^)rj, is given in (5.32). With this, we find that the attenuation coefficient is given 

by 



Here we recall that (Aa;^(t))rp and {X^{t))rj. are given in (5.32) and (5.33). 
Now, for times short compared with the Ohmic decay time we find 

kT^2 

aT(t) ^ exp{ ^ r2 cP}, t < m/C. (5.43) 

This is exactly the form of the attenuation coefficient for a free particle [7] , which for very 
short times is of the form aT{t) = exp{— i^/rj}, where the decoherence time is 

r. = (5.44) 



in which v — ^JkT/m is the mean thermal velocity. 



VI. CONCLUDING REMARKS 



The system we are discussing is that of an oscillator coupled to a linear passive heat 
bath, with a microscopic Hamiltonian of the form (2.9). The long time equilibrium state of 
this system is entangled, in the sense that the normal modes correspond to coupled motion 
of the oscillator and the bath. The Hu-Paz-Zhang equation is an exact master equation 
describing how this entangled equilibrium state arises from an initial state in which the 
bath and the oscillator are not coupled. We should perhaps emphasize that this assumption 
of a decoupled initial state is common to all derivations of a master equation, going back 
at least to the work of Wangness and Bloch [19], who phrased it as an assumption that at 
any instant of time the system is approximately decoupled. Indeed, such an assumption is 
essential for the introduction of the notion of partial trace, i.e., the trace over states of the 
uncoupled bath as in (4.1), key to the existence of any master equation. Now, our exact 
solution has allowed us to see more clearly how serious is that assumption. In particular, 
we have seen that within this assumption an exact solution leads to meaningful results only 
in the high temperature limit. Here we hasten to add that this remark does not apply to 
the many successful applications of master equations in weak coupling approximation. 

For the most part, previous discussions have been made under the restriction that the 
initial state of the particle is a pure state, effectively at zero temperature, while the bath is 
at a high temperature T. There has even been an approximate experimental realization of 
such a state [20] . However, we would argue that such a state is unphysical in the sense that 
it does not correspond to what is envisioned when one speaks of a system at temperature T. 
Rather, the initial time dependence is then dominated by the "warming up" of the particle, 
which occurs on a time scale of order the decay time m/C,. On the other hand, as we have 
shown, the restriction to such a state is not necessary, one can, within the assumption of an 
uncoupled initial state, choose the particle state to be at the same temperature as the bath. 
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In order to describe a state of the system that is entangled at all times, including the 
initial time, it is necessary to abandon master equation methods. Some time ago, a more 
general method applicable to such systems was described by Ford and Lewis [21]. In their 
method, a system in equilibrium is put into an initial state (e.g., a wave-packet state) by a 
measurement and then at a later time is sampled by a second measurement. This method of 
successive measurements has recently been applied to obtain exact results for the problems 
of wave packet spreading and decoherence [6]. For the wave packet spreading one finds in 
place of (5.22) or (5.31) the result 

where the variance is now given by 

= + (6.2) 

An equivalent result for wave packet spreading in the Ohmic case has been obtained by 
Hakim and Ambegaokar [17], who used functional integration methods. For the decoherence 
problem, one obtains in place of (5.38) or (5.42) the result 

s(t) 

Note that these results are finite at any temperature and apply to an arbitrary heat bath 
and for all times. At short times, the results (5.34) and (5.43) are in agreement with these 
exact results. 



APPENDIX A: SINGLE RELAXATION TIME MODEL 

Here we consider the single relaxation time model for the case of a free particle {K = 0). 
This model corresponds to a memory function of the form 

pi{t) = ^e-'N{t), (Al) 
r 

where 6 is the Heaviside function. Note that in the limit r — > this becomes the Ohmic 
memory function (5.1). With this form of the memory function and with K — 0, the response 
function (2.7) takes the form 

a{z) 



—mz^ — iz , 

I — IZT 



—mz{z + i'y+jiz + i7_) ' 
where we have introduced 



^ + ^(7+ + 7-) ^A2) 



4Cr 
m 



7± = — (A3) 
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Note that in the Ohmic hmit 7+ — > — > 00 and 7_ — > C/m and we recover the form (5.2) 
of the response function. 

With this form of the response function, the Green function (2.6) can be written in the 
form 



7^(1 -e-^-*)- 7^1 



e 



G{t) = ^ \ ' ' ^, (A4) 

^7- 7+ (7+ -7-) 



If we form the quantities (2.19) with this Green function, we find 

(7+ - 7-) [(7+ + l-)e^-^ - 7+] 



2r{t) = 7- - 



(7+ + 7_) (eT+* - eT-*) + 7+ - 7_ ' 



^2.^N 7-7+(7+ - 7-) /^gx 

(7+ + 7_)(eT+*-eT-*)+7+-7_' ^ ^ 

In the Ohmic limit, it is clear that for any finite time, 2r(t) 7_ ^ C/""^ ^''^^ Q'^{t) 0. 

On the other hand 2r(0) = and ^2^(0) = 7-7+ — ^ 00. For t of order r we see that 

in this Ohmic limit 

^yt) ^ ^e-'l^ ^ ^5(t). (A6) 
mr m 

This justifies the assertion made in (5.4). 

Next, we consider the fluctuating position operator, X(t), for the single relaxation time 
model. With the Green function (A4) it is a simple matter to verify the relation 

G(-t') - G{t - f) = -mG{t)G{-t') - i:i^^^^G(-t'). (A7) 

7+7- 

Putting this in (5.5) we find in place of (5.7), 

Tfl'T 

X{t) = x,{t) - x,{0) - mG{t)x,{0) - —[1 - mG{t)]xs{0). (A8) 
Forming the mean square, we can write 

X\t)) = s{t) - mG{t)s{t) + -m''G\t)s{Q) 



c 



2 2 

^ [l-mG{t)Ys^'^\Q). (A9) 



2^ 

The interest here is in the zero temperature limit. With the response function given by (A2) 
and with 7 = 0, the expression (5.6) for the mean square displacement can be put in the 
form 



sit) = ^ '^'''-rr''^\ (Alo) 



2fl 7^/(7_t)-7^ /(7+t) 
ttC 7+ - 7- 

where I{x) is given by (5.14). We see now that s(0) is finite, but the fourth derivative s(^)(0) 
is logarithmically divergent. Indeed the divergence is the same as in the Ohmic case, with 
the same long time form (5.17) 
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